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Abstract Several recent advances in turbulent dynamo theory are reviewed. High res- 
00 ■ olution simulations of small-scale and large-scale dynamo action in periodic domains 

O^ ' are compared with each other and contrasted with similar results at low magnetic 

Prandtl numbers. It is argued that all the different cases show similarities at interme- 
diate length scales. On the other hand, in the presence of helicity of the turbulence, 
power develops on large scales, which is not present in non-helical small-scale turbulent 
dynamos. At small length scales, differences occur in connection with the dissipation 
i-i^ I cutoff scales associated with the respective value of the magnetic Prandtl number. 

»— ^1 These differences are found to be independent of whether or not there is large-scale 

O ' dynamo action. However, large-scale dynamos in homogeneous systems are shown to 

f-j . suffer from resistive slow-down even at intermediate length scales. The results from 

tj^ ' simulations are connected to mean field theory and its applications. Recent work on 

helicity fluxes to alleviate large-scale dynamo quenching, shear dynamos, nonlocal ef- 
fects and magnetic structures from strong density stratification are highlighted. Several 
insights which arise from analytic considerations of small-scale dynamos are discussed. 
> 
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1 Introduction 

Dynamos convert kinetic energy into magnetic energy. In the astrophysical context, 
one always means by a dynamo a self-excited device, except that the conductivity is 
uniformly distributed and not limited to conducting wires in an insulating container, as 
in laboratory dynamos. Since the beginning of the space age, it is well understood that 
such a devi ce can work, at least in principle. The discovery of the first rigorously proven 
example bv lHerzenberd ( 1958t ) was significant, even though the particular configuration 



studied by him was not of immediate astrophysical interest. The Herzenberg configu- 
ration consists of a uniformly conducting solid medium (e.g., copper) with two or more 
conducting rotors in electrical contact with the rest. Dynamo action is possible when 
the rotors spin faster than a certain critical value that depends on the angle between 
the rotors. For angles between 90 and 180 degrees the solutions are non-oscillatory, 
while for angles between an d 90 degrees there are o scillatory solutions that were 
discovered only more recently I Brandenburg et al.lll998l 'l. 



The Herzenberg dynamo belongs to the class of "slow" dynamos, for which the 
dynamo growth rate is maximum for intermediate values of the conductivity, but it 
returns to zero in the limit of perfect conductivity. Another example of a slow dynamo is 
the Roberts flow, which consists of a steady two-dimensional flow pattern, u = u{x, y), 
but all three f l ow co mponents are non-vanishing. This dynamo was first studied by 
[Robert j (|l970l . Il972l l. It provides an important benchmark of a large-scale dynamo, 
where magnetic field on scales larger than the scale of motions are produced, by a 
mechanism called the a effect. Here, a refers to the name of the coefficient in the 
relation between mean magnetic fi eld B and me an electromotive force f in a turbulent 
flow. The basic i dea goes back tolParkeiJ ( 1955|) who proposed a relation of the form 



£ — olB. Later, ISteenbeck et al.l (196a) computed a tensorial relation of the form 
£i — ctijBj for rotating stratifled turbulence. It was clear that higher derivatives of 
the magnetic field would also enter, so a more general expression is 

£i = aij o Bj +ri,jk° Bj^k, (1) 

where the isotropic part of the tensor ry^Tj, jyjjj, — r/teijfei corresponds to turbulent 
magnetic, diffusion where -qt is the turbulent magnetic diffusivity. The circles indicate 
convolution over time and space, which is however commonly replaced by a multipli- 
cation in the limit in which th e integral kernels are w ell approximated by 5 functions 
in space and time. The work of ISteenbeck et al.l ( 19661 ) marked the beginning of mean- 
field electrodynamics, which is still the leading theory to explain the amplification of 
magnetic fields on length scales that are larger than the scale of the energy-carrying 
turbulent eddies. Such systems are therefore also referred to as large-scale dynamos. 
Essential here is that the trace of the a tensor is non- vanishing. This typically requires 
helicity of streamlines (at low conductivity) or of vortex lines (at high conductivity) . 
In this review we discuss our current knowledge of dynamos covering both large- 
scale and small-scale turbulent dynamos. In practice, all dynamos of astrophysical 
relevance tend to be "fast" and thus have a finite growth rate also in the limit of perfect 
conductivity. There are examples of predetermined fiows for whic h fast dynamo action 
is ind icated by numerical simulations, for ex ample the ABC fiow ( Gallowav and FrischI 



1986l l and the Gallo wav and Proctor! (1992) fiow, but the convergence of the a effect 



with increasing magnet ic Reynolds number can be quite slow or non-existent (see, e.g., 
ICourvoisier et al.ll2006r ). By contrast, when relaxing the assumption of a predetermined 
kinematic flow pattern and adopting a turbulent fiow field, the a effect appears to be 



converged for magnetic Reynolds numbers exceeding a cr itical value of t he order of 
unity up to values of about 200 probed in the simulations ( Sur et al.ll2008l ). 



We begin by discussing examples of numerical realizations of turbulent dynamos 
and then turn to some astrophysical examples. Large-scale dynamos produce fields 
that are well characterized by suitable averages. A theory for describing the evolution 
of such averaged fields is mean-field theory that is obtained by averaging the governing 
equations, most notably the induction equation. Large-scale dynamos can then also be 
referred to as mean-field dynamos (MFD). Under certain conditions, making suitable 
assumptions about the spatial variation of the a effect, solutions of the mean-field in- 
duction equation have been used to characterize the large-scale magnetic fields seen in 
astrophysical bodies like the Sun and some spiral galaxies. Mean-field theory is also ap- 
plied to the momentum equation. This leads to a number of effects inclu ding turbulent 
visco sity, the A effect (responsible driving differential solar ro tation; see lRiidigerlll980l . 
[l989|), the anisotropic kinetic alp ha effect ([Fris ch et al.l ll987l ). as well as the negative 



effective magnetic pressure effect ( Kleeorin et al.lll989t ). In this paper, we also discuss 



important effects resulting from the mean-field momentum equation, namely the spon- 
taneous formation of magnetic flux concentrations in a strongly stratified layer. Such 
results may be relevant to explaining the formation of active regions in the Sun. 

There is now increased interest in what is known as small-scale or fluctuation 
dynamos. These ar e dynamos t hat ca n work alre ady under fully isotropic condi t ions and 
were anticipated bv lBatchelon ( 1950!) and others (Biormann and Schliiter 195ll : lElsasserl 
[l956), but the now accepted theory was provided by Kazantsev ( f967 ). Fluctuation 
(or small-scale) dynamos are important in cosmic objects because they are generic to 
any random flow of a sufficiently conducting plasma. Furthermore, the growth rate is 
fast and random fields can grow on the eddy turnover timescale of the smallest eddies, 
which are typically much shorter than the age of the system. This is particularly true of 
galaxy clusters, for which small-scale dynamos are crucial to explaining the observed 
magnet ic fields, because the con ditions for large-scale dynamo action are probably 
absent ( Subramanian et al.ll2006l ). 



Small-scale dynamos are nowadays also invoked to describe the small-scale magnetic 
field at the solar surface. However, in many contexts both small-scale and large-scale 
dynamos go together. And then it is not clear whether one can (or should) distinguish 
between a small-scale field from a small-scale dynamo from that associated with the 
fluctuations that are inherent to any large-scale dynamo and that can be caused by 
tangling and amplification of the large-scale field. 



2 Numerical realizations of turbulent dynamos 

Next, we discuss some numerical realizations of turbulent dynamos. By "turbulence" we 
mean here fiows that are solutions to the Navier-Stokes equations that are irregular in 
space and time, subject to energy injection at large length scales and energy dissipation 
at short length scales. The ratio between forcing scale to dissipation scale is quantified 
by the Reynolds number. Throughout this paper, we adopt the following definition for 
this number: 

Re = Urms/fkf, (2) 

where Urms — {u ) ' is the root mean square (rms) value of the turbulent velocity, 
v is the kinematic viscosity, and fcf is the forcing wavenumber, i.e., the wavenumber 



of energy injection. A similar definition is adopted for the magnetic Reynolds number, 
Rsm = Urms/rikf, where -q is the magnetic diffusivity. We also define the magnetic 
Prandtl number, Prjvf — v/rj = Rejv,//Re. 

At small Reynolds numbers, Re <^ 1, the flow is determined by viscosity and has 
qualitatively different properties from flows with Re ^ 1. The break point is close 
to unity, so, for example, a flow with Re = 5 would be called turbulent (or at least 
mildly turbulent), because it begins to show certain asymptotic scaling properties that 
are also found for fully turbulent flows. The most famous aspect of fluid turbulence 
is the Kolmogorov E}^{k) ~ k~ ' energy spectrum, which is normalized such that 
J E}^ dfc — Urms/2. Density is omitted in this definition, which requires that density 
fluctuations are unimportant and the flow is nearly incompressible. In this review we 
also discuss stratifled flows in which density varies strongly due to gravity. In those 
cases the incompressibility condition V ■ n = has to be replaced by V • pu = 0, which 
is valid as long as the flows are slow compared with the sound speed and the typical 
scales of variation smaller than a scale height. However, in the following we do not 
make such assumptions and consider fully compressible flows. 



2.1 Small-scale and large-scale dynamos 

We begin by discussing first the difference between small-scale and large-scale dynamos. 
A small-scale dynamo is one that generates magnetic field at scales much smaller than 
that of the energy-carrying eddies, while a large-scale dynamo generates field at scales 
larger than that of the energy-carrying eddies. In may practical situations, the difference 
is of somewhat academic interest, because the small-scale dynamo is always excited 
when the magnetic Reynolds number is large, which is the case in many astrophysical 
situations. Furthermore, the conditions for the excitation of large-scale dynamo action 
are met in many situations of interest. However, in the case of isotropic turbulence 
such a distinction can be made by considering helical and non-helical turbulence. In 
both cases the system can be homogeneous and it makes sense to compute spectra of 
magnetic and kinetic energy, E-M{k) and E^{k), respectively; see Figure[T]for Prjv/ = 1. 
The early evolution of such a dynamo is quite similar: both dynamos h a ve a k ' power 
spectrum at small scales. Suc h a scaling was predicted bv iKazantsevI 1 1967h (see also 



iKulsrud and Anderson! (119921 )) for a single scale non-helical flow which was 5-correlated 
in time, but seems to be obtained in the simulations more generally. A t late times, 



howe ver, helical turbulence allows the development of an inverse cascade IjFrisch et al.l 
I1975I ) on a longer resistive timescale (see below). 

We see that, at least in the saturated state, large-scale dynamos produce and sus- 
tain magnetic fields at scales larger than the energy injection scale, while small-scale 
dynamos produce and sustain magnetic fields at scales smaller than the energy injection 
scale. The lack of similar behavior in the linear regime could be interpreted as evidence 
that th e underlying mechan ism for producing large-scale fields must be nonlinear in 
nature ([Herault et al.ll2011f ). Alternatively, one could interpret the resulting dynamo 
as a combination of large-scale and small-scale dynamo action, where the latter has 
a larger growth rate such that in the kinematic regime the field is dominated by the 
small-scale field, although the large-scale dynamo does still operate in the background. 
As a consequence it gets the chance to dominate only when the small-scale dynamo 
has already saturated. Even in this case, and if the Re^v/ is large, we will argue that 
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Fig. 1 Kinetic and magnetic energy spectra in a turbulence simulation without net helicity 
(up) and with net helicity (bottom) for a magnetic Prandtl number of unity and a mesh size is 
512"^ meshpoints. Notice the pronounced peak of Ej^ik) at fc = fci in the case with helicity. The 
energy input wavenumbers are k{ = l.Sfci in the non-helical case (upper panel, Rejv/ = 600, 
PfM = 1) and fcf/fci = 4 in the helical c ase ( lower plot, Rcm = 450, Pr M = 1). Adapted from 
iBrandenbure and SubramanianI ll2005af) and iBrandenburg et al.1 | |2008D . 



one needs to additionally dissipate small-scale magnetic helicity before the largest scale 
field can appear. 

Both large-scale and small-scale dynamos show that the spectral magnetic energy 
exceeds the kinetic spectral energy in the beginning of the inertial range. This dif- 
ference is slightly weaker for large-scale dynamos, but this might be an artifact of 
the Reynolds number still not being large enough. Such a difference us ed to be com- 
pletely absent at Reynolds numbers previously reported; see, for example IBrandenburg 
(|200ll ). The slight super-equipartition was quite e vident w hen simulations at a reso- 
lution of 1024 meshpoints became first available ([ Haugen ct al. 2003|). a l though this 
feature can already be seen in earlier simulations (jMeneguzzi et al.lll98ll : iKida et al.l 
Il99ll ). This spectral excess of magnetic fields is expected to diminish as one pro- 
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Fig. 2 Magnetic and kinetic energy spectra for runs with 512'* meshpoints and hyperviscosity 
with hyperresistivity (solid line) and Sniagorinsky viscosity with hyperresistivity (red, dashed 
line). Note the mutual approach of kinetic an d magnetic energy sp ectra before entering the 
dissipative subrange. Adapted from lHaugen and Brandenburg ||2006| ). 



ceeds furthe r down the inertial range. Such behavior was indeed seen in simulations of 
Haueen and Brandenburg (2006 ) when using a combination of Smagorinsky viscosity 
for the velocity field and hyperresistivity for the magnetic field; see Figure [21 This 
implies that the two spectra cannot be parallel to each other at intermediate length 
scales, and that the slope of Eyi{k) must be slightly steeper than that of E]^{k). This 
difference in the two slopes at intermediate wavenumbers is now asso ciated with the 
observ ed differences in the spectral exponents in the solar wind; see iBoldvrev et al.l 
(201l|), who find steeper spectra for the magnetic field than the velocity field both 
from simulations and solar wind observations. 

f n the spectra of Figure [T] we see a remarkable difference between small-scale and 
large-scale dynamos. In particular, large-scale dynamos are capable of producing a peak 
of magnetic energy at the smallest possible wavenumber, while small-scale dynamos 
do not. On the other hand, we have stated earlier that small-scale and large-scale 
dynamos are difficult to distinguish in the early stage. In Figure [3] we show the critical 
values of Rejv/ (= Re) for dynamo action both for small-scale (non-helical) and large- 
scale (helical) dynamos. Note that for Rej\,/ > 35 the growth rate attains a Re ' 
scaling both for small-scale and large-scale dynamos. The same growth rates a re also 
obtained for dynamos driven by convection; see Fig. 15 of iKaovla et alJ l|2008l ). This 
Re ' scaling of the growth rate of the rms magnetic field, impUes that the growth 
rate is not controlled by the turnover time of the energy- carrying eddies, but of eddies 
at the dissipation scale ( Schekochihin et al.ll2002l . l2004al ). At low values of Rejv/, only 
large-scale dynamo action remains possible. Its excitation condition is determined by 
the requirement that a certain dynamo number exceeds a critical value. This usually 
translates to t he condition that th e degree of scale separation is large enough; see 
equation (5) of lBrandenburei ( 2003 ). 



An important issue for astrop hysical applications is how coherent are the fields 
generated by small-scale dynamos ( Subramanian et al.ll2006l : ICho and Rvull2009l ). On 
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Fig. 3 Dependence of dynamo growt h rates of the rms m agnetic field on Rejy/ for helical and 
non- helical turbulence. Adapted from iBrandenburd ||2009| ). 



the basis of simulations done with large Prjv^, but small Re. lSchekochihin et al.l ( 2004b|) 
argued that the field saturates with a folded structure, where the field reverses at the 
folds such th at power concentrates on resistive scales. The simulations of lHaugen et al.l 
(2003, I2OOJ) with Rcm = Re ^ 1, found the magnetic correlation lengths ~ l/6th the 
velocity correlation l ength, but much larg er than the resistive scale. This seems consis- 
tent with the simple ISubramanianI (1999) model for nonlinear saturation of small-scale 
dynamos. What happens at large Re and large Ptm, which is representative of galactic 
and cluster plasmas, is not easy to capture in simulations. Th e Pim — 50, Re — 80 
simulation described in Branden burg and SubramanianI (l2005a), showed evidence for 
strong field regions with folds, but equally, regions with strong fields and no folds, il- 
lustrating that such structures need not be volume fill ing. Moreove r, the 'spontaneous 
stochasticity', that applies to highly turbulent fiows ( Evinkll201ll ). suggests that the 
dynamics of small-scale dynamos could be quite different in turbulent compared to 
laminar high-Pr/yf systems. Furthermore, in galaxy clusters the viscosity may be set 
by plasma effects (jSchekochihin et aL 2005al'l. It could also be highly anisotropic owing 
to the presence of magnetic fields i Parrish et al.ll2012f l. In addition, heat condition is 
also very a nisotropic, giving rise to magnetothermal and heat fl ux-driven buoyancy in- 
stabilities ( Parrish and Stoneil2005l : |Parrish and Quataertll2008l '). Further work on these 
aspects is desirable, using both semi-analytical ideas and high resolution simulations. 



2.2 Low magnetic Prandtl number 



In many astrophysical settings, like solar, stellar, or accretion disk plasmas, the mag- 
netic Prandtl number is rather low (~ 10~ or less). This becomes numerically hard 
to handle, especially if the magnetic Reynolds number should still be large enough 
to support dynamo action. In Figure |4] we compare non-helical and helical runs with 
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Fig. 4 Kinetic and magnetic energy spectra for non-helical (top) and helical (bottom) turbu- 
lence at low magnetic Prandtl numbers of Prjvf = 0.02 and 0.01, respectively. Here, Rejvf = 230 
with Re = 11,500 in the non-helical case, adapted from Brandenburg (2011a), and Rej^f = 23 
with Re = 2300 in the helical case adapted from iBrandenburg (2009 ) . In both cases the 
resolution is 512^ mesh points. 



magnetic Prandtl numbers of 0.02 and 0.01, respectively. In the former case the mag- 
netic Reynolds number Rejv/ is 230, which is weakl y supercritical ; the c ritical value 
of Rejv/ for small values of Pi m is ^m.crit ~ 150 ( Brandenburg l2011al l , compared 
with /?m,crit ~ 35 for Prjvf = 1 (Haugen et al. 20 04|). These value s agree with those of 
earlier work (|Schekochihin et al. 2005b . 2007 : Iskakov et al.lbOOTi '). So, for Rbm = 230 



and Prj\,/ = 0.02, we have Re = Re^/PrM = 11,500. Normally, this would require a 
numerical resolution of about 10, 000 meshpoints, but it turns out that at low values 
of Pr^,/, most of the energy is dissipated resistively, leaving thus very little kinetic 
energy to be cascaded, terminating therefore the kinetic energy cascade earlier than at 
PrAf = 1. 
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Fig. 5 Visualizations of Bx/Baq on the periphery of the domain at six times during the late 
saturation stage of the dynamo when a large-scale field is gradually building up. The small- 
scale field has reached its final value after t/r ^ 100 turnover times. Here, _Beq = V^OPO ^^rms 
is the equipartition field strength where kinetic and magnetic energy densities are comparable, 
and po is the mean density. Note that the maximum field strength is about twice Beq. 



At intermediate length scales, kinetic and magnetic energy spectra are close to each 
other. The magnetic energy spectrum no longer exceeds the kinetic energy spectrum, 
as was found for Pr^j/ = 1; see Figure [T] Again, this might be a consequence of still 
insufficiently large Reynolds numbers and limited resolution. In fact, it is plausible 
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Fig. 6 Saturation of the small-scale magnetic energy density and continued increase of the 
large-scale magnetic energy density. Here, kf/ki = 15. 



that, in the limit of large fluid and magnetic Reynolds numbers, kinetic and magnetic 
energy spectra coincide, even if Prj;/ is small. And only at much smaller scales, the 
magnetic energy spectrum turns into a dissipative subrange, and goes below the ki- 
netic energ y spectrum, due to s tronger Ohmic dissipation. The slopes of the k~ ' 
spectrum of Golitsvnl (Il960ll and Moffatti (Il961f) and the scale-invariant fc~ spectrum 



( Ruzmaikin and Shukurovll982l : iKleeorin and Rogachevskiilll994l : iKleeorin et al.ll 19961 1 
are shown for comparison. However, the kinetic energy spectrum tends to be steeper 



than k 



-5/3 



and is closer to k 



-1.9 



and k' 



-7/3 



in the non-helical and helical cases. 



2.3 Inverse transfer and a effect 



The dynamics of a large-scale dynamo is most dramatic when the scale separation 
ratio is large, i.e., fcf/fci ^ 1. In Figure [T] it was only 4, but now we consider a case 
where kf/ki — 15. In Figure [5] we show visualizations of one component of the field 
for different times. Evidently, a large-scale field develops that varies in the y direction. 
This particular large-scale field is best described by xz averages. In Figure [6] we show 
the evolution of the mean energy density of this large-scale magnetic field, {B ), and 
compare it with that of the small-scale field, {b ) = {B ) — {B ). Note that the small- 
scale field reaches its final saturation value during the time span considered, while the 
large-scale field has not yet saturated and is expected to do so on the diffusive time 
scale of the box such that trjki = 0(1). It is also interesting to note that the large- 
scale field starts becoming noticeable only when trjki ~ 0.02, or when (t/t^) ~ 4.5, 
where ij = {rikf)~ is the resistive timescale at the forcing scale. In other words one 
needs several resistive (diffusive) times at the forcing scale, before the large-scale field 
can grow. We will interpret this result below in terms of the resistive dissipation of 
small-scale magnetic helicity which alleviates a quenching. 
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Fig. 7 Magnetic energy spectra Eyii^'), at earlier (top) and later (bottom) times. The scale 
separation ratio is kf/ki = 15. The range of time t is given in units of the turnover time, 
T = 1/nrmsfcf. At small wavenumbers, the Eyi{k) spectrum is proportional to fc^, while to the 
right of kf/ki = 15 there is a short range with a k~^'^ spectrum. 



The evolution of the magnetic energy spectrum for this case is shown in Figure [3 
where we see several stages during the early phase of the dynamo, and in a separate 
panel the later saturation behavior. Clearly, a large-scale field is already present for 
t/r > 100, but the field is then still fairly isotropic and therefore not very pronounced 
in visualizations shown in Figure [S] Spectra of magnetic energy and rescaled magnetic 
helicity are shown in Figure |8] for the saturated state. Here, magnetic helicity spectra 
Hy[(k) are normalized such that J Hyi dfc = (A-B), where i? = V x A is the magnetic 
field expressed in terms of its vector potential. Note that at early times, the magnetic 
field shows the Kazantsev k ' slope (t/r = 10) in the range 7 < k/ki < 25. However, 
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Fig. 8 Spectra of magnetic energy, Eyi{k), and rescaled magnetic helicity, ■iikHyi{k)/2^Q 
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Fig. 9 Left panel: time dependence of the peak wavenumber for scale separation ratios of 
15 (dashed) and 30 (solid lines) at ReM of 12, 27, and around 57 (increasing line thickness). 
Right panel: Rejv/ dependence of the cascade speed for scale separation ratios of 15 (open 

symbols) and 30 filled symb ols. The straight lines give the Re^jj (dotted) and Re^j^ (dashed) 
dependences. Adapted from iBrandenburtj ll2011bf) . 



already at times t/r = 20 and 30 one sees a small hump at fcf/2, which is significant 
in view of an interpretation of these results in terms of a so-called a dynamo, which 
will be discussed in Section [2.41 

The temporal increase of the typical scale of the large-scale field can be determined 
by monitoring the quantity 



lm{t) = k~\t)= f k-^EM{k)dk/ f EM{k)dk. 



(3) 



In Figure [5] we plot k^^ (t) for different values Rejvf- Here, we normalize with respect 
to the minimal wavenumber fci = 27r/L in the domain of size L, so we plot fci/fem- 
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There is a limited scaling range that allows us to determine the increase of £m(t) for 
different values of Rejv/ . We measure the speed at with the bump travels toward larger 
scales by the quantity Qtiav = d£m/dt. In the right-hand panel of Figure [^ we plot 
the resulting values of atrav as a function of Rejvf . The results are compatible with a 

-1/2 

resistively limited inc rease whose speed diminish es like Re^.y^ . This behavior was first 



seen in simulations of [Brandenburg et al.1 ( 2002 ) and then more convincingly at larger 



resolution in simulations of iBrandenburd ( 2011a ). Such a behavior further reinforces 



our earlier remark that the large-scale field can only grow on the slow resistive timescale 
in periodic boxes. 



2.4 Connection with mean-field theory 



There exists a close analogy between the inverse transfer described above and mean- 
field dynamo theory in that both are able to predict the occurrence of large-scale fields 
with similar properties and excitation conditions. In mean-field theory one splits the 
velocity field U and magnetic field B into the sum of a mean, large-scale components 
(C/ and B) and a turbulent, stochastic components (n and b); that is C7 = C7-|-it and 
B — B + b. One then solves the averaged induction equation, 

^^Vx{UxB + £-vf^oJ), (4) 



where £ = u x b is the mean electromotive force that we discussed already in connec- 
tion with Equation ([l|. Under the assumption of isotropy and sufficient scale separation 
in space and time, we have just £ = aB — rjtjJ.oJ, where a and rjt are a pseudo-scalar 
and a scalar respectively. For the case when there is no mean flow, a stability analysis 
gives the dispersion relation for the growth rate A as 

\ = ak- {rit + rj)k , (5) 

and the eigenfunctions are force-free solutions with rit^oJ — aB, which are plane 
polarized waves, just like in Figure [S] where the large-scale field can be approximated 
by B oc (sinfc2/,0, cos fcj/), ignoring here an arbitrary phase shift in the y direction. 
The dynamo is excited when Ca = Q/(??t + ??)fc > 1, where Ca is the relevant dynamo 
number in this context. The fasted growing mode occurs at wavenumber k = fcmax = 
a/[2(r;t + v)]- Furthermore, using estimates for the high-conductivity limit, 

a = Ok = "^''"('^ ■ ''*) 8'iid Vt ~ ■|''"('" )) (6) 

we find that Ca = effcf/fci ( Brandenburg et al.ll2002l ). where tf < 1 is the fractional 



helicity, and rj <^ rjt has been assumed. We can now return to our discussion in con- 
nection with Figure [71 where we notice that at early times the field growth occurs at 
wavenumber A:f/2. This is indeed the value expected from our simple estimate, since 
for fully helical turbulence, we expect {ijJ-u) — k[{u ), and thus fcmax ~ a/2rit = fef/2. 
It is of interest at this point to comment on the validity of mean-field dynamo 
concepts. In condensed matter physics for example mean-field theory is generally valid 
when applied to systems where fluctuations are assumed small. In high Re a/ turbulent 
systems on the other hand, the fluctuations grow more rapidly than the mean-field, due 
to small scale dynamo action. Thus, even in the kinematic stage when Lorentz forces 
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are small, one needs a closure theory to calculate the mean-field turbulent coefficients 
like a and rjt- Traditionally these coefficients have been derived by still making a quasi- 
linear approximation (strictly valid for small fluctu ating fields), which is also known 
as the first order smoothing approximation ( FOSA) l Mofi^attlll978l : lKrause and Radlerl 



ll98Ci : [Brandenburg and Subramanianll2005al ). This is sometimes also referred to as the 
second order correlation approximation. Some improvements to this can be made by 
adopting closure approximations like the minimal tau approximation whereby triple 
correlations involving fluctuating flelds are taken to provide a damping term pro- 



portional to the quadratic correlations ( Blackman and Fieldll2002l : iRadler et al.ll2003l : 



[Brand enburg and Subramanian 2005a). There are also a few cas es, like S correlated 
flows tKazantsev 1967; Zeldovich ct al. 1983) or renovating flows ( Dittrich et al.lll984l : 
[Gilbert and Bavlvll 19921 ) which provide analytically solvable models, where the form of 
a and rft given by Equation (O is recovered. 

In this context, direct simulations as discussed above, which can be interpreted in 
terms of mean-field concepts, lend some validity to the theory. This applies also to 
the interpretation of results in the nonlinear regime to be discussed below. Moreover, 
when the a and rjt have been measured directly in simulations of isotropic turbulence, 
one gets results remarkably close to the estimates of FOSA given in Equation © 



( Sur et al.ll2008l ). This suggests that the strong magnetic field fiuctuations produced 



by small-scale dynamo action do not contribute a systematic large-scale component to 
the mean emf £ — u x b, correlated with the mean field. They do make the mean- 
field coefficients noisy. However, the fact that we can still use the mean-field concept 
in understanding the results of direct simulations implies that this noise does not 
have a crucial effect, perhaps after the small-scale dynamo has nearly saturated. The 
saturation of the dynamo will be discussed further in Section [ 



2.5 Shear dynamos 

Remarkably, not all large-scale dynamos require an a effect. In fact, large- scale dynamo 
action has been seen in simulations with just shear and no helicity; see iBrandenburd 
( 2005al ) for simulations using a shear profile motivated by that of th e Sun. An obvious 



cand i date a t the time was the so-called shear-current effect ( Rogachevskii and KleeorinI 
2003I . I2OO4I ). which requires riijUij > 0, where rjij is the part of the magnetic diffu- 
sivity tensor that multiplies J such that £i — ... — rjijjj, and Uij is the velocity 
shear. Already early calculations using the test-field method showed that the relevant 
component of rjij has the wron g sign (IBrandenburg 2005b| ). This confirmed the re- 



sults of quasilinear calculations (IRadler and Stepano- 



20061 : iRiidiger and Kitchatinovl 



I2OO6I : ISridhar and Singhll20ld : ISingh and SridhaiJl201ll ). Moreover in the large ReM 
limit using FOSA, but for arbitrarily strong shear, the co rresponding cross coupling 
implied by th e shear current effect was shown to be absent I Sridhar and SubramanianI 
I2009all2009bl ). The issue of how the mean- field grows in nonhelical turbulence in the 
presence of shear remained open in view of other possible contenders. 

One possibility is the incoherent a-shear dynamo that lives fr om the combination of 

shear and fiuctuations of the a effect and was originally invoked bv lVishniac and Brandenburg 
l| 19971 ) to explain the magnetic field of unstratificd shearing box simulations of accretion 
disc turbulence (Hawlcv ct al. 1995). This mechanism has received considerable atten- 
tion i n the follow ing years (Sokolov 19 9?!: ISilant'evlbOOOl: iFedotov et aklbood : IFroctoiJ 



20071 : iKleeorin and Rogachevskiii i2008l : ISur and SubramanianI |2009| ) . Meanwhile 



ev- 
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idence for the existence of shear dynamos in simple shearing box simulations wa s 
mounting JYousef et all l2008all2008bl : iBrandenburg et alJ boosi : iHerault et aD boilh . 



Although the underlying mechanism may have appeared to be a new one, there is now 
quantitative evidence that this can be explained by an incoherent g-shear dy namo 

This 



quantitative eviaence tnat tnis can pe expiamea Dy an mconerent g-snear gy n 
( Brandenburg et al.ll20oi : iHeinemann et al.l[20nl : iMitra and Brandenburg||2012[ i. ' 



is remarkable given the unconventional nature of the approach whereby one uses mean- 
field theory over two spatial directions and considers the fluctuations that remain in 
time and the third coordinate direction as physically meaningful. 



2.6 g-quenching 

A fully satisfactory theory for the magnetic feedback on the g effect does not ex- 
ist. What we do know is that for strong mean fields B, this a becomes a tensor 
of the form Oij = ai{B)5ij — a2{B)BiBj, where B = B/B are unit vectors, and 
B — \B\ is the modulus. However, if this tensor is applied to the B field, we have 
otijBj = («! — a2)Bi, which suggests that a is just like a scalar. We also know 
that at least part of the quenching acts in such a way that the total field (small- 
scale and large-scale) obeys the magnetic helicity ev olution equation. Th is was de- 
rived some time ago in a certain approximation by IPouQuet et al.l ( 119761) and was 
then applied to derive an equation for the quenching of a (IKleeorin and Ruzm aikin 



1%2; Gruzinov and Diamond 19 93: iKleeorin et al.l I2OO0I: 



Field and Blackman 



Blackman and Brandenburdl2002l : ISubramanianll2002l : IBrandenburg and Subramanian 
2005a| ). 



The crucial starting point is the realization of lPouquet et al.l ( 197(a ) that under the 



infiuence of Lorentz forces, the g effect has an additional component, gjvj = hrj ■ b/pQ, 
where j ■ bis the current helicity associated with the small-s cale field and a — UK +aM 
is the sum of kinetic and magnetic a effects. Interestingly. IPouquet et al.l ( 19761 ) also 



showed that rit does not get renormalized under the same approximation. Under locally 
isotropic conditions, in the Coulomb gauge, j ■ b can be approximated by k^a- b/fiQ, 
where a ■ b = hf is the magnetic helicity of the small-scale fields. In order to write an 
evolution equation for the magnetic helicity density one can fix a gauge for the vector 
potentia l. One could also work in terms of the evolution equation for the current 
helicity (|Subramanian and Brandenburg! 12004 ). Perhaps more elegant is to write this 
evolution in terms of a gauge invariant magnetic helicity density, defined as the density 
of corr elated links of b, and which is most closely related to the h{ in the Coulomb 
gauge ( Subramanian and Brandenburg||2006l ). The evolution equation for hf is 



^ = -2e-B-27^kfhf-V-Ff, (7) 

where Ff is the magnetic helicity fiux of the small-scale field. This equation shows that 
the g effect produces magnetic helicity at a rate —2£ ■ B = —2a^^^B , where g,od ~ 
a — r]tkm is the reduced a effect and km ~ fJ-aJ ■ B/B is the effective wavenumber 
of the mean field. In a supercritical dynamo, the sign of gj-od agrees with that of a 
(the rjt term is subdominant). Then starting with a specific sign for the kinetic gx and 
zero magnetic gjvi, this produces gjvi of opposite sign, which quenches the total g and 
the dynamo progressively with increasing field strength. In the absence of a magnetic 
helicity fiux, this process happens on a resistive time scale, which is what is seen in 
Figure [6l where final saturation is not even remotely in sight. We recall that a rapid 
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Fig. 10 Magnetic energy spectra Ej^(k) of positively (red) and negatively (blue) polarized 
parts at earlier times. Note the preferred build-up of the large-scale field at E^(k{/2). relative 
to the slower growth of E^ (fcf ) • Slopes proportional to fc-^ for _E^ and fc* for E~ are shown. 



evolution of the energy of the me an field up to fci /fcf times th e equipartition value is 
expected on theoretical grounds i Blackman and Brandenburg| |2002). In practice, this 
is hard to verify because at early times the mean field has not yet reached the scale 
of the system and modes of different orientation are still competing. Nevertheless, by 
splitting the magnetic field into its positively and negatively polarized contributions, 



fields (IBrandenburg et all I2OO2I : IBi 



it_ is possible to separate large-sca le and small-scale 

Brandenburg and Subramanian 2005al : iBrandenburd 

l2011a ). In Figure [10] we clearly see a faster build-up of the large-scale field through 
E7,{k) compared with the small-scale field through EZ,{k). As we have argued before, 
this build-up of large-scale fields is still resistively slow, but it is important to realize 
that the demonstrated existence of large-scale fields in the kinematic stage provides 
support for the usefulness of the mean-field approach. 

The final saturation value, in periodic box simulations, can be estimated simply by 
noting that in the absence of magnetic helicity fiuxes, the total current helicity must 
vanish in the steady saturated case, i.e., J ■ B + j ■ b — 0. This is because the total 
current helicity drives the change of magnetic helicity, which in steady state must be 
zero. Such a state can be obtained in a nontrivial manner with helical forcing when 
current helicity has opposite signs at large and small length scales. For example, if the 
kinetic helicity at small scales is positive (that is qk is negative), then the generated 
qm is positive. In that case the current helicity of small-scale fields is also positive with 
Moi • b = fcf b^ and hence (iqJ ■ B — —kiB is negative. (This implies that fcm ~ — fci-) 
Assuming furthermore equipartition betw een kinetic and ma gnetic energies at small 
scales, i.e., b^ ~ Mopi*^ 



Boq, we obtain ( Brandenburgj|200ll ') 



SV-Be'q 



fcf/fcl 



(8) 
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in the final state. We recall that in the run shown in Figure [S] the scale separation 
ratio is kf/ki = 15, so it is understandable tha t there is not yet any sign of saturation 
of the large-scale field; see iBrandenburg ( 200ll ) for early results on the resistively slow 



saturation. 

We do not expect resistively slow saturation beh avior to occur in star s or gala xies , 
because the V-.Ff term can usually not be neglected jBlacknian and Field L2000a 20 00bl : 



Li 

b|2 



iKleeorin ct al. 2000:iBra ndenburg and Subranianiari2005bl : Shukurov et al. 2006). These 
results were obtained by solving the mean-field equations. Subsequent simulations have 
shown that tur bulent diffusive fluxes exist that would constitute a s ufficient contribu- 
tion to 'V ■ F( ( Mitra et al.l l20ld : iHubbard and Brandenburg |2010|) , especially when 



Rejvf is larger than a critical value around 10 . This flux term can then dominate 



over the 2rij ■ b term ( Candelaresi et al.ll2011l ). Under certain considerations it is pos- 



sible that a flux of the form £ x A from the mean electromotive force contributes to 
the flux term, but by solving an evolution equation for the magnetic helicity density 
of the total (small- and large-scale) field, consideration of this term can be avoided. 
This approach is also suited to deal with fluxes ass ociated with gauges that can intro - 



duce artificial fiuxes in shearing environments; see lHubbard and Brandenburg ( 20 111 ) 
They find no evidence for a fiux resul ting from shear that were previously argued to 
be important (1 Vishniac and Chq 2001 ) . Using simulations with anisotropic non-helical 



forcing in the presence of shear, IShapovalov and Vishniad (|201l|) argue that large-scale 
dynamos might even live entirely due to helicity fluxes, although the exact origin of 
this flux remains to be clarified. Another natural contribution to the flux term is just 
advection of both small-scale and large-scale fields, along with the assoc iated mag- 
netic helicity ( Shukurov et al.ll2006l : ISubramanian and Brandenburg l2006h . This will 



naturally arise from coronal mass ejections in the solar context, or supernovae driven 
outfiows in galajcies (IBlackman and BrandenburdboOa : IShukurov et ahlboOa : ISur et al-l 
l2007l : IWarnecke et al.ll2oTil )! 



Large-scale dynamos therefore seem to need helicity fluxes to work efficiently. This 
conclusion can be understood more physically as follows. As the large-scale mean field 
grows, the turbulent emf £ is transferring helicity between the small- and large-scale 
fields. The large-scale helicity is in the links of the mean poloidal and toroidal fields 
of the astrophysical system like the Sun or the Galaxy, while the small-scale helicity 
is in what can be described as 'twist helicity' (or simply twist) of the small-scale field, 
produced by helical motions. Lorentz forces associated with the 'twisted' small-scale 
field would like to untwist the fleld. This would lead to an effective magnetic qm effect 
which opposes the kinetic ax produced by the helical turbulence. The cancellation of 
the total a effect can lead to catastrophic quenching of the dynamo. This quenching 
can be avoided if there is some way of dissipating the twist (which is slow in high-ReM 
system) or transferring the twists in the small-scale field out of the region of dynamo 
action, or cancelling it between two hemispheres. That is, if there are helicity fiuxes 
out of the system or between different parts of the system, the large-scale field can 
grow to observable strengths. 



2.7 Application of mean-field theory to galaxies 

The mean-field theory described above has been applied extensively to understand 
magnetic fields of disk galaxies. The mean-field dynamo equations allow substantial 
simplification provided a suitable parameterization of turbulent transport coefficients 
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is chosen. Of course, this parameterization presumes a suitable closure for nonlinear 
effects to arrive at a closed set of nonlinear mean-field dynamo equations. Such an 
approach does not necessarily imply a deep understanding of the physical processes 
involved in the magnetic field evolution. However, it appears to be sufficient for prag- 
matic modeling of magnetic field configurations in particular galaxies to be compared 
with the observational data of polarized radio emission. 

A first example of such parameterizat ion and further drastic simplification of mean- 
field dynamo equations was presented bv lParkeil (195^). Early simplified models for the 
ga lactic dynamo, which a llow analytic or quasi-analytic investigations, ca n be found 



Ruzmaik in et al. (198 q), and more re cent reviews bv iBeck et al.l 1 19961 ) . IShukurovl 



(200J), and Kulsrud an d Zweibell ( 20081 ) . Next, one can suggest, as the most prag- 



matic contender, simple mean-field model fo r galactic dynamo the so-called no-z model 
( Subramanian and Mestell 119931 : lMosslll995l ). The idea of this model is to present az- 



imuthal and radial components of the mean galactic magnetic field by their quantities 
at the galactic equator and average the mean-field equations with respect to the co- 
ordinate z perpendicular to the galactic plane. The third component of the magnetic 
field can be reconstructed from the condition divB = 0. This approach is an obvious 
oversimplification, but it appears still adeq uate for modeling m agnetic field evolution 
including the helicity fluxes discussed above (|Sur et al.ll2007l : IChamandv et al. 2012'). It 
also allows one to model magnetic conflgurations for design studies of new gen erations 
of radio telescopes such as the Square Kilometer Array; cf. IMoss et al.l 1 20121 ). 



2.8 Application to mean-field dynamos of stars and the Sun 

In comparison with galactic dynamos, the st atus of the theory of s olar a nd stellar 
dynamos is less satisfactory. Early models by ISteenbeck and Krausd ( 1969|) provided 



numerical solutions to the full two-dimensional axisymmetric mean-field equations un- 
der the assumption of an assumed profile of the internal angular velocity fi(r, 6'), whose 
radial gradient, dfi/dr, was negative, and variations in colatitude 6 were ignored in 
most of their models. This yielded cyclic magnetic fields with equatorward migra- 
tion under the assumption that the a effect is predominantly positive in the northern 
hemisphere. Subsequently, helioseismology delivered detailed contours of 0(r, 9), which 
excluded the previously assumed H profiles and suggested that dQjdr > at low lati- 
tudes where strong magnetic flux belts are observed to propagate equatorward during 

the course of the 11 year sunspot cycle. 

Various solutions have been offere d to this dilemma l|Parkerll 19871). The most promi- 
nent one is the flux transport dynamo i Choudhuri et al.ll995l : lDikDati and Charbonneaul 



[l999|), whereby meridional circulation causes the dynamo wave to turn around in the 
opposite direction. This model operates under the assumption that the turbulent mag- 
netic diffusivity operating on the toroidal fleld is much lower than the value suggested 
by mixing length theory, and the a effect is assumed to work only at the surface. 
The other proposal is that the d ynamo wave obeys equatorward migration owing to a 



narrow near-surface sh ear layer ( Brandenbur g 2005ai '). where dQjdr is indeed strongly 
negative; see Fig. 4 of iBenevolenskava et al.l ( 19991 ). This proposal still lacks detailed 



modeling. In view of the shortcomings in the treatment of mean-field dynamo theory 
(e.g., our ignorance concerning nonlinearity discussed in Section [2.61 or the neglect of 
finite scale separation discussed below in Section [2. lOp . the ground for speculation re- 
mains fertile. Magnetic helicity fluxes in interface and flux transport dynamos have 
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already been looked at ( Chatteriee et alJuOlOh . but finite scale separation effects are 



neglected. Indeed, the Sun is strongly stratified with its scale height changing rapidly 
with depth, making it hard to imagine that simple-minded approaches that ignore this 
can be meaningful at all. A cornerstone for the proposal that solar activity is a shallow 
phenomenon may lie in the success of explaining the formation of active regions and 
perhaps even sunspots a s a result of spontaneous form ation of flux concentrations by 
convective fiux collapse [Kitchati nov and Mazuril200Cl ) or the negative effective mag- 
netic pressure instability; see Section [2.111 below. 



2.9 Magnetic helicity in the solar dynamo exterior 

There is now explicit evidence for the presence of magnetic helicity in the exterior 
of dynamos. In particular, it has been possible to detect magnetic helicity of oppo- 
site signs at small and large length scales. This has been possible through measure- 
ments of magnetic helicit y spectra in the solar wind away from the equatorial plane 
( Brandenburg et al.ll201ll ). This data pointed for the first time to a reversal of magnetic 



helicity density between interior and exterior of the dynamo. Such reversals have now 
been confirmed in numerical s imulations of dy n amos coupled to an exterior. There are 
first of all the simulations of IWarnecke et al.l l|201ll ') showing coronal mass ejections 
from a turbulent dynamo in a wedge of a spherical shell, where the reversal occurred 
some distance away from the dynamo. Next, there are related simulations in Cartesian 
geom etry where a revers al can be f ound immediately above the surface; see figure 12 
of Warnecke and Brandenburg ( 2010 '). Finally, there are earlier me an-field simulations 
showing such a reversal as well; see the lower panel of figure 7 of [Brandenburg et al.l 

Hooi). 

The occurrence of such a reversal is now well understood in terms of the magnetic 
helicity equation for the small-scale field shown in Equation (0. In the dynamo interior, 
the Q effect dominates over turbulent magnetic diffusion and produces magnetic helicity 
of a sign opposite to that of a. In the northern hemisphere, a is positive, so this 
produces negative h{. Turbulent magnetic diffusivity opposes this effect, but it is still 
subdominant and can therefore not change the sign of h{. This is different in the solar 
wind where the a becomes subdominant compared with turbulent diffusion, which is 



the reason the sign of /if is now different IjBrandenburg et al.ll201ll 



2.10 Scale separation in space and time 



In connection with Equation ^ we noted that the relation between £ and B does, 
in general, involve a convolution in space and time. This becomes important when the 
variations of B occur on time and length scales comparable with those of the turbu- 
lence, whose turnover time is l/urms and its typical scale is ^ = fcr . The properties 
of the integral kernel are often determined in Fourier space, in which case a useful 
approximation of it is &ij{k, cj) — a^- K{k,Lj) and 'i)^jj.{k,u>) = r; .!K{k,uj), where 

K{k,ij)^- . ^ ,„^„ . (9) 

^ ' l-liUT + k^l^ ^ ' 

Such an integral kernel has rece ntly been obtained with the test-fi eld method ap- 
plied to passive scalar diffusion ( Rheinhardt and Brandenburg! 12012 ), and in limit- 



ing cases (either with u = or with A; = 0) for a and r]t\ see [Brandenburg et al.l 
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(12003) and iHubbard and Brandenburg (|2009l ') for details and applications to spatial 
and temporal nonlocalities, respectively. The test-field method allows one to deter- 
mine the turbulent transport coefficients by solving an extra set of equations that 
describe the evolution of the fluctuating magnetic field for each test field, which is a 
predetermined mean field. Under some conditions it is also necessary to solve corre- 
spon ding evolution equations for velocity perturbations i Rheinhardt and Brandenburg 



[2010). The combined presence of sp atio-temporal nonlocality was first considered by 
(jRhcinhardt and BrandenburgbOlZ ). who proposed Equation ([5]) and used it to refor- 
mulate Equation ([l]) as a simple differential equation of the form 

^ + "I ^ ^'£) ^* = "'?^^ + '^S^^.'^- (10) 

Such an equation is quite easy to implement and represents an improvement in terms 
of physical realism. This representation avoids not only the problem of causality associ- 
ated with the infinite spee d of signal propagation in the absence of the time derivative 



( Brandenburg et al.l 12004 1 . but it also prevents the development of artificially small 



scales in the mean field. 

The application of this new technique is still in its infancy, and it needs to be seen to 
what extent spatio-temporal nonlocality can substantially alter the nature of the solu- 
tions. As an example we note that a finite r has been fou nd to lower the critical dynamo 
number for oscillatory solutions by a factor of about 2 ( Rheinhardt and Brandenburg 



I2OI2I ). In the context of disk galaxies, such non- locality in time can also lead to phase 
shifts between the spiral forcing of the dynamo by matter arms, and the resulting 
magnetic spirals, as seen in the galaxy NGC6946 ( Chamandv et al.ll2012l 'l. 



2.11 Magnetic structures resulting from strong density stratification 

Finally, let us discuss a mean-field effect that occurs under the condition of strong den- 



f inally, let us discuss a mean-neld eltect tnat occurs under tne co ndition ot strong den- 
sity stratification. It has been theoret i cally anticipated long ago (IKleeorin et al.l 19931. 



1996 : lKleeorin and Rogachevskiil994l : lKleeorin et al.ll989l . ll99(]l : lRogachevskii and KleeorinI 



20071 ) ■ but only more recently was it also seen in numerical simulations of the mean 



field equations ( Brandenburg et al.l 2010p a nd then in direct numerical simulations 



( Brandenburg et al.ll201ll : iKemel et al.ll2012l l 



The essence of this effect is the suppression of turbulent intensity by a mean mag- 
netic field. This means that the effective pressure caused by B is not just B /2fiQ, but 

there must also be an additional cont ribution from the suppression of the turbulence, 

which leads to a negative contribution (J Kleeorin et al.[l989| . l l99d : iKleeorin and Rogachevskiil 
I1994I : iKleeorin et al]ll996l : lRogachevsku and KleeorinI [20071 1. This modifies the nature 
of the magnetic buoyancy instability in such a way that magnetic structures become 
heavier than their surr oundings and sink. This has been demonstrated using both 
mean- field simulations ([Brandenburg et al.ll20ld ') as well as direct numerical simula- 



tions ( Brandenburg et al.ll201ir ). As an example we show in Figure ITT] a snapshot from 



a direct simulation where a weak (-Bo/^oqO = 0.01) magnetic field is imposed in the y 
direction. The density stratification is isothermal, so the density scale height is constant 
in the direction of gravity (the negative z direction). The total density contrast from 
bottom to top is about 540. The instability grows at a rate which scales with rjtQk , 
where 7)to = Urms/Skf is an estimate for the turbulent ma gnetic diffusivit y, which is 
well reproduced by simulations using the test-field method ( Sur et al.ll2008l ). 
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Fig. 11 Visualizations of By for a simulation with strong density stratification and a weak 
imposed magnetic field (0, BoiO) with Bo/BcqO = 0.01, Rejvf = 18 and Ptm = 0.5. Time is 
given in units of turbulent— diffusive times. Note the gradual emergence of a large-scale mag- 
netic flux con centration wh i ch th en sinks as a result of negative effective magnetic buoyancy. 
Adapted from lKemel et al.l ||2012|) . 



The study of this negative effective magnetic pressure is still very much in progress. 
In particular, it has not yet been studied how this negative effective magnetic pressure 
instability interacts with the mean-field dynamo. It is envisaged that this instability 
might produce local magnetic field enhancements in the surface layers that resem- 
ble active regions. Such regions are long lived (~ 1/2 year). Traditionally such long 
timescales have not been associated with the surface layers. However, the time scale 
of the negative effective magnetic pressure instability is the turbulent-diffusive time, 
'''td = ('?to'^i)~ 1 where 77to ~ '"rms/3fcf is the estimated turbulent magnetic diffusivity. 
This time can be much longer than the local turnover time, rto = (wrmsfcf)" , which 
is about (1 day at the bottom of the near-surface shear layer at 40 Mm depth. The 



3(fcf/fci 



which can be around 300 for a scale 



ratio of these time scales is Ttj/rto 
separation ratio of just 10. 

In the rest of this paper we focus on small-scale magnetic fields that occur over 
a range of different astrophysical settings. They are believed to be important in un- 
derstanding small-scale magnetic fields in the surface layers of the Sun, although that 
part of the field might also be a consequence of tangling the large-scale magnetic field. 
Another possible astrophysical application of small-scale dynamos might be the clus- 
ters of galaxies, because in the absence of rotation it is difficult to motivate any form 
of large-scale dynamo action. 
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3 Analytical approaches to small-scale turbulent dynamos 

Direct numerical simulations are now the most straightforward way to understand tur- 
bulent dynamos. However, more traditional analytical method provide a useful support 
for them. Analytic methods have provided particular insights both into large-scale and 
small-scale dynamos. We discuss some specific analytic considerations of small-scale 
dynamos further below. 



3.1 Correlation tensor and small-scale dynamo 

A natural approach here is to introduce the second-order correlation tensor of the 
magnetic field Bi{x,t) as 

Bijix,y,ti,t2) = {B,{x,ti)Bj{y,t2)), (11) 

taken at two spatial points x and y and at two instants ti and t2- Here (...) denotes 
averaging over an ensemble of turbulent velocity field fluctuations which can be de- 
scribed by the velocity fleld correlation tensor Vij{x,y,ti,t2) constructed in the same 
way as Bij. 

Fo r a particular model of turbulence (short-correl a ted r andom flow), iKazantsevI 
(|l967l ). and simultaneously iKraichnan and NaearaiarJ (|l967l ') for a slightly different 
model, obtained a gov erning equation known now as the Kazantsev equation. In par- 
ticular |Kazantsev| (1965) assumed Vij (a;, i/,ii,f2) = Vij{x,y)S{ti—t2), and derived an 



evolution equation for the magnetic fleld correlation tensor calculated at two simulta- 
neous instants Bij{x,y,t) — Bij{x,y,t,t), with ti — t2 ~ t. This reads 

where Lj^ji^^ is a second-order differential operator with coefficients depending on V^j , 
its spatial derivatives and coefficient of magnetic diffusion rj. In some sense, the Kazant- 
sev equation is similar to the famous Steenbeck-Krause-Radler a effect equation in 
mean-field electrodynamics. In practice however the latter equation provided much 
more astronomically fruitful results than the first one. The reason is presumably two- 
fold. First of all, the Kazantsev equation requires more algebra for its solution than 
the mean-field equations. We only take here some n ew points isolated rece ntly in this 
bulky algebra and refer to a detailed review given bv lZeldovich et al.l (|l99d l. The other 
reason is that the physical interpretation of solutions of Eq. p2[) is more delicate than 
that for the mean-field equation. This is the main issue presented in following sections. 
Original insight which lead Kazantsev to Eq. (|12p required some mathematics 
from quantum field theory that can hardly be considered 'user friendly' for a per- 
son with an ordinary MHD background. A more familiar approach, which is how- 
ever rather bulky, can be foun d elsewhere ( Zeldovich et al.l Il990l : ISubramanianI 119971 : 



[Brandenburg and Subramanian 2005a ) . The particular form of the Kazantsev equation 
(|12[) is associated with a specific model of turbulence assumed, in particular the as- 
sumption of 5-correlated (in time) velocity fiuctuations. This model is quite restrictive 
and does not allow one to represent adequately some basic properties of the Kol- 
mogorov cascade. Various attempts to derive th is equations for more realistic models 
were undertaken (see, e.g.. iKleeorin et al.ll2002l ). For example, incorporating a finite 
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correlation time results in quite complicated albeit beautiful mathematics and gives 
rise to integral equations. As far as it is known however the results are more or less 
the same if just applied to Eq. (|12[l with Vij taken from a suitable model of turbulence 
ignoring the fact that this very model is incompatible with the derivation of Eq. (|12p . 
A reasonable way to simplify Eq. H12[) to a useful level is to consider statistically 
homogeneous, isotropic and mirror-symmetric turbulence and look for solutions having 
the same symmetry. This means that we exploit a velocity correlation tensor of the 
form 

Vy = A{r)rirj + B{r)5,j (13) 

{r = y — x) and look for the magnetic field correlation tensor in a similar form. 
The incompressibility condition means that the functions A and B depend on a single 
function, say, F{r) while the solenoidality condition of the magnetic field means that 
Bij depends on a single function, say, W{r,t). Then Eq. (|12p can be reduced to a single 
second-order ordinary differential equation for a single function which depends on W 
with coefB c ients depending on F. In fact it is the equation which was obtained by 
iKazantsevI (jlOGT]). The algebra here remains quite bulky and we avoid to prese nt it 
here in detail; see, for example , the detailed discussions in lZeldovich et al.l (1990|) and 
[Brandenburg and Subramania n (2005a). 

There is no problem to solve the Kazantsev equation in the homogeneous and 
isotropic case for a particular choice of F numerically or by analytical approximations. 
In fact the Kazantsev equation can be reformulated as a Schrodinger type equation 
for a particle with variable mass m{r) in a potential U{r) (both of which depend on 
F{r)), whose bound states correspond to exponentially growing W{r,t). It is possible to 
develop a WKB-like method for an approximate solution in the limit of large magnetic 
Reynolds numbers Re^ (~ ?;~ ). 

A general result following from these solutions can be summarized as follows. For 
a sufficiently large Rejv/ , the magnetic field correlation tensor (i.e. W) grows exponen- 
tially with a growth rate 72 which is determined by l/v where I is the turbulence corre- 
lation time and v is its rms velocity. The critical magnetic Reynolds number is of order 
-Rm.crit ~ 10 , with _Rtii crit = 26 in the most simple example with F = exp—{r/l) . 
For this illustrative example, function W which determines the magnetic field correla- 
tion properties has quite a complicated form which contains spatial scales from / up to 
ZRe-i/^ 

A plausible scenario for the nonlinear saturation of the growth of the magnetic field 
correl ation tensor governing the Kazantsev equation was suggested by ISubramaniarJ 
(|l999l ). The main idea here is that nonlinearity results in an effective suppression of 
Rcjvf up to i?ni,crit- It is quite straightforward to consider Kazantsev equation for ho- 
mogeneous, isotropic and mirror-asymmetric turbulence and to c ombine concepts of 
the second-order correlation tensor with the mean-field approach ( Subramanianlll999l : 



iGabov and Sokolofll l2004') . Interestingly, the Kazantsev equation in the presence of ki- 
netic helicity, can be reformulated into a tunnelling type quantum mechanical problem, 
where by the bound states of the small scale dynamo can 'tunnel' to develop long range 
correlations ( Subramanianll 19991 : iBrandenburg and Subramanianll200d : iBoldvrev et al.l 



|2005). It is also possible to solve the Kazantsev e quation for locally homogeneous and 
isotropic turbulence in a finite body of size L ^ I ( Maslova et al.ll 19871 : iBelvanin et al.l 
[l993i ). 
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3.2 Random magnetic fields in cosmology 

The description of a random magnetic field by a second - order correlation function is 
widely exploited in cosmology. To be specific, IChernirJ (1966|) considered at first a 



cosmological model with a random magnetic field and discussed cosmological evolution 
of its magnetic energy. We have to stress however that the concepts of statistical 
hydrodynamics, which are a mathematical basis for the Kazantsev approach, need 
some modification to be applicable for curved spaces of General Relativity. The point 
is that Eq. (|lip considers the product of two vectors applied at different spatial points 
while such option is absent in Riemannian geometry. We have to consider a geodesic 
line connecting points x and y, transport a field from the second point to the first one 
and consider product of the two fields applied at one spatial point. This recipe presumes 
that a geodesic lines which leaves a point x does no longer cross other geodesic lines 
which leave this point. In other words this means that geodesic lines do not contain 
conjugated points and the space-time has no gravitational lenses. If they do exist 
we have to elaborate somehow the concept of the correlation tensor and no general 
recommendations are suggested until now. 

Then statistical homogeneity, isotropy and mirror symmetry in a curved spatial 
section of a Friedmann cosmological model reads 

e*^' = C{r)n'n^ + D{r)g^\ (14) 

i.e. we have to distinguish upper and lower indices and use tangent vector n' of the 
geodesic line connecting x and y instead of the vector r which do not exist in a curved 
space. The other point is that one deals with a curved space formulation of solenoidal- 
ity condition to reduce C and D with one one function. Inspired by earlier work of 
iGarcia de Andradd ( 2010l ). iRubashnv and S okolofi (2010) performed a corresponding 



analysis to show that for the Lobachevsky space with negative curvature, 

C^-R — ^F', D = F + R — 2.i?'. (15) 

where R is the curvature radius of the spatial section and r is distance between x and 
y. It is instructive to compare this representation with that for Euclidean space 

C = ~^F', D = F+'-F'. (16) 

Because th -rr has a finite limit at r/R — >■ oo, correlations decay slightly slower for the 
Lobachevsky space rather in the Euclidean space (of course F is the same in both cases). 
More significantly, however, the volume of a sphere with radius r grows exponentially 
in the Lobachevsky space and not as a power law like in Euclidean space. It means 
that F has to decay much faster in Lobachevsky space than in Euclidean space to get 
convergence of various spatial means which are based on the correlation tensor. 
For a spherical space (closed cosmological model) one obtains 

C = -R^1^F', D = F + R^-^^F' . (17) 

This representation gives CijiR) = 0, i.e., C vanishes if y is just the opposite point to 
X. This looks reasonable because n* is not determined uniquely for the opposite points, 
tg -^ diverges for r = -rR/2 so the finiteness of correlations means that F'{^) — 0. 
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This condition is specific for spherical geometry and has no direct analogues for the Eu- 
clidean one. Remarkably, a closed universe admits another homogeneous and isotropic 
topological structure of spatial section, namely an elliptical space which is twice smaller 
than the spherical one. Instead of the above condition we find here F" {■kR/2) = 0. 
Moreover, an elliptical space does not admit orientation, so one cannot distinguish be- 
tween left- and right-hand coordinate system there. It means that quantities such as 
helicities, a effect and other pseudo-scalar quantities cannot be introduced there. Such 
topological constrains on the magnetic field correlation properties look rather strange. 
Fortunately they do not affect substantially physically interesting conclusions because 
the correlation length / is usually much smaller than the curvature radius R. The prob- 
lem however is that both quantities as well as the horizon radius vary strongly during 
the course of cosmological evolution that is given I being negligible in the present day 
cosmological scale might be very large in scales of the Early Universe. A more severe 
problem arises if we are goi ng to consider a homogene ous and isotropic ensemble of 
random gravitational waves ( Ivanova and Sokoloffll2008l '). which are often discussed in 
cosmological context. 



3.3 Higher statistical moments and intermittency 

To get more detailed information concerning a dynamo-generated small-scale magnetic 
field, it is useful to consider higher statistical moments which are introduced as en- 
semble averages of a product of p magnetic field vectors (p is the order of statistical 
moment). Following the Kazantsev approach one can obtain the governing equations 
for these quantit ies and demonstrate that the moments grow provided Re^y,/ is high 
enough (see, e.g.. lKleeorin et al.ll2002l ). Of course, the algebra becomes more bulky as 
m increases. The problem is that the higher moments grow faster than the lower ones 
in the sense that 

72/2 < 74/4 < 76/6 . . . (18) 

Of course, this fact can be supported by a direct calculation. However, it is much more 
instr uctive to demonstrate the phenomenon at a qualitative level fe.g. Molch anov et al.l 
1 19881 ). Let us consider a flow with a memory time r so the magnetic fleld B{nT) at 
instant nr can be considered to be developed from the initial fleld B(Q) which is 
affected by n independent random transport operators Ti 

B(nr) = f„f„„i...TiB(0). (19) 



As a matter of fact, progressive growth of higher statistical moments in Eq. (|18|) does 
not depend critically on the fine structure of operators Tj so we can illustrate the 
phenomenon by considering the simplest operators Tj, i.e. just independent random 
numbers Ti. 

For the sake of definiteness, let InTj have Gaussian distribution with zero mean 
and standard deviation a. Then 

T = TnTn-l...Ti (20) 

is a log-normal random quantity and In T has zero mean and standard deviation ^/na — 
ay/t/T. A straightforward calculation shows that 

(TP) = exp(aVi/2r), (21) 
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so the normalized growth rates of the moments 7p/p ~ cr p/2t grow Unearly with the 
degree p of statistical moment. The other point is that the growth is determined not 
by a typical value of Tj which is of order a, but by strong deviations which are of order 
aWtp/r. The probability density of a Gaussian quantity to achieve the level aWtp/r 
is of order exp(— ip/2r) so the size A^ of a statistical ensemble should be as large as 
N*{t) = exp(tp/2r) to include such rare events. 

Note that the above analysis presumes that the statistical ensemble is infinitely 
large. If the ensemble is only large but finite its size A'' should exceed a critical value 
N*{t) which grows in time exponentially. For large t we obtain A^ < N* (t) and the 
above estimates becomes inapplicable. 

If we consider a medium of independent cells of size I renovating after a memory 
time r, then the critical size of the system which allows one to recognize the growth 
governed by Eq. H'ill) is given by 

L*{t) = IN^^^ = ;exp(fp/6r). (22) 

This means that the behavior of the statistical moment of the order p is determined 
by very rare cells and Eq. (|22|) gives an estimate for the distance to the nearest cell 
which determines the moment at a given point. The phenomenon of a random field 
whose properties are determined by rare and remote events is known as intermittency. 
The wording comes from medicine and means a state when a person is near death, 
but his/her heart still works from time to time. These rare events of the heart activity 
determine the fact that the person is still alive. 

Note that if we calculated a PDF of T based on a limited sample with A'^ < A'^* (i) it 
is practically impossible to recognize the existence of the above mentioned rare events 
which do not contribute to the PDF calculated. Of course, the importance of the result 
depends on how large t should be to make the intermittency recognizable and how 
large the corresponding value A'^* is in comparison with A'^ typical of celestial bodies. 

It is natural to address this point based on a simple physical example rather than 
just a product of random operator s. A simple example of this kind accessible for simple 
numerics has been suggested bv IZeldovichI (196J) in a cosmological context. Let a 



remote object have a (small) angular size 9, let x be the distance to the object and 
y = 6x its linear size. y(x) is known in the Riemannian geometry as Jacobi field and 
is governed by the so-called Jacobi equation 

y" + Ky = 0, (23) 

where a prime mean s the derivative taken with respect to x and K is the spatial 
(sectional) curvature. IZeldovichI (196J) recognized the importance of density and then 



curvature fluctuations on the evolution of y along the line of sight. In other words, we 
consider K as a. random, say, Gaussian quantity. 

It is quite easy to simulate many independe nt solutions of Eq. I|23[) and dete rmine 



experimentally how large N*{t) is. It appears ( Artvushkova and SokolofB l2005l l that 
one needs N ^ 5 x 10 . . . 10 to recognize the difference between 72/2 and 74/4 for 
i ~ 10 . Of course, a simulation of 10 independent 3D cells for a hundred turnover 
times becomes prohibitive given the purely computational problems, in addition to 
the problems associated with the further data processing of the results. On the other 
hand, the number of independent turbulent cells in a typical galaxy is of order of 
10 so the effects of intermittency can contribute in mean quantities of interest for 
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galactic dynamos. In practice however the difTerence between 72/2 and 74/4 is for the 
intermittent fields, as far as it is known, not very large. 

The other point here is that the growing magnetic field becomes sooner or later dy- 
namically important. Of course, it is important to know which happens first, whether 
the field becomes dynamically important first or whether the size of the domain be- 
comes too small to contain the intermittent structure? Unfortunately, the abilities of 
simple models like Eq. (|23|) to reproduce the stage of nonlinear dynamo saturation are 
limited. A numerical experiment with a simple model shows that the number A'^ re- 
quired to reproduce the behavior of higher statistical moments declines strongly when 
the solution becomes dynamically important such that exponential growth of statistical 
moments saturates. 

Yet another point which can be recognized from the experiences with simple models 
is that e nsemble averagin g is not the only option to describe the behavior of the growing 
solution. iLamburt et al.l (,2003) demonstrated that the quantity 

-^^ -^ 7 > (for i ^ 00) (24) 

for almost all realizations (with probability 1, or "almost sure" in the wording of 
probability theory) where no averaging is taken at all. Quantities such as 7 are known 
as self- averaging quantities and 7 in particular is known as Lyapunov exponent. The 
phenomenon becomes clearer if one introduces a 2D row- vector {y, y') and rewrites 
Eq. pS)) as a vectorial equation 

z' = zA, (25) 

where 4 is a random matrix process with vanishing trace (i.e. Tr A = 0). Then the evo- 
lution of z from an initial condition zq can be represented as a product of independent 
random unimodular matrices Bi = exp(^ir), detB^ = 1, where Ai is a realization at 
a given interval of renovation of the random matrix process A. 

The product of independent random mat rices is quite we l l inve stigated in prob- 
ability theory (so-called Furstenberg theory). IZeldovich et al.l (|l984l ) stressed the im- 
portance of this theory for small-scale dynamos (here the elements of A are dvi/dxj). 
iMolchanov et al.l (|l984f ) argue that magnetic field generated by a small-scale dynamo 
grows such that 

lnlB|(a;,t) , , , 

— ' — '^^^ ^ 7 > for t^ CO, (26) 

where 7 is a positive constant. A numerical experiment for the Jacobi equation supports 
this interpretation and shows that t should be of the order of a hundred memory times 



to get this behavior (jArtvushkova and Sokolofill2005l) . T he PDF of the dynamo-excited 
magnetic field is investigated bv lChertkov et al.l l|l999h . 



Unfortunately, this approach of simple models cannot mimic the a effect. The point 
is that one cannot produce a pseudoscalar quantity a based on correlations of Aij and 
Amn and a Levi-Civita tensor epqr- 

Generally speaking, the analytical results discussed above show that detailed direct 
numerical simulations of small-scale dynamos at the kinematic stage can be nontriv- 
ial to interpret due to the rapidly growing intermittency. Fortunately, the magnetic 
field becomes dynamically important quite rapidly (at least at small scales), so the dy- 
namo becomes nonlinear and mathematical distincti ons between pro perties of various 
statistical moments become less and less important. lArtvushkova an d Sokolofi ( 2 0051 ) 
investigated several simple models how catastrophic intermittency typical of the kine- 
matic stage gradually evaporate when nonlinear effects become important. Presumably, 
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something like that happens for the much more complicated full small-scale dynamo 
equations. 



3.4 Small-scale magnetic field and shell models 

The Kazantsev model is developed for the kinematic or weakly nonlinear stage of the 
dynamo and its ability to describe the strongly nonlinear stage of a dynamo is obviously 
quite limited. However, it remains useful to have a simple model of strongly nonlinear 
dynamos in terms of ordinary (instead of partial) differential equations. An option of 
this kind is given by so-called shell models of MHD-turbulence. 

The starting point of the shell model approach is to note that for random fields and 
flows, one can hardly reproduce in a numerical simulation the actual realization of the 
random field which are obtained in a given celestial body. In practice we are interested 
in some spectral properties of the field of interest. If we can get such properties without 
solving the full set of equations, we would be happy with such a result. 

Shell models are designed to describe the cascade process over a large range of 
scales (wavenumbers) by a chain of variables Un(t), bn(t), each of them characterizing 
velocity or magnetic field oscillations with wavenumbers k in the range from fc„ = fco-^" 
to fc„_|_i i.e., a shell of wavenumbers. The parameter A characterizes the ratio of two 
adjacent scales (the width of the shell) and usually A < 2. The model includes a 
corresponding set of ordinary differential equations, which should reproduce the basic 
properties of the equation of motion. In particular, the model has to reproduce the type 
of nonlinearity of the primitive equations and to retain the same integrals of motion in 
the dissipationless limit. Let us note that shell models can possess positively defined 
integrals of motion (energy, enstrophy in two-dimensional turbulence, and the square 
of magnetic potential in 2D MHD-turbulence) , as well as quadratic integrals with an 
arbitrary sign (the integrals of this kind are usually called 'helicities'). The signs of the 
helicities are defined by the balance between the contributions of odd and even shells 
to corresponding quantity. 

The shell m odels were sug gested by Kolmog orov's school to describe the spectral 
energy transfer ( Gledzerll 19731 : Desnianskii and Novikovlll974l ). After numerous refine- 



ments they became an effective tool for description of the spectral properties of the 
small-scale turbulence (see for review Bohr et al. 199o). The shell models for M HD tur- 
bulen ce we re introduced by [F rick (1984), Gloagucn ct al. (1985), Brandenb urg et al.l 
(|l996r ). and lFrick and Sokolofa (11998! ). This approach reveals many intrinsic features of 
small-scale dynamo action in fully developed turbulence of co nducting fluids (for review 
see lBif eralc 2003). In particular, the shell model suggested bv lFrick and Sokolofll (|l998l ) 
gives a fast growth of small-scale magnetic flelds (on the timescale l/v) and its satura- 
tion at the equipartition level as well as non-Gaussian (similar to lognormal) PDF for 
small-scale magnetic field in the saturated state. In some cases shell models give a hint 
concerning dynamo actio n in the parametric dom ain inaccessi ble to dir ect numerical 
simulations. In particular. IStepanov and Pluniani (2006) and iFrick et al.l (|2006 'l argue, 
based on simulations of MHD-shell models, that the critical magnetic Reynolds num- 
ber for small-scale dynamo action remains moderate -Rm.crit ~ 80 in the case of low 
Prandtl numbers. This is also confirmed by the simulations discussed in Section 12.21 
General speaking, shell models seem to provide an effective way to investigate small- 
scale dynamos. 
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It may be possible to combine shell models as a tool to describe small-scale vari- 
ables in a dynamo with grid or spectral methods for large-scale variables in mean-field 
equations. However on l y the first steps in th is direction have been made until now 
( Frick et al.|[20o3 . bood : JNigro and Veltri|[20lil . 



3.5 Dynamical chaos and small-scale dynamo 

One more point of comparison between numerical and analytical approaches to the 
small-scale dynamos is as follows. The Kazantsev model and its related investigations 
consider turbulence as a truly random field and apply concepts of probability theory 
in full extent. In contrast, direct numerical simulations and shell models consider tur- 
bulence as a chaotic behavior of solutions for deterministic equations of motion. It is 
quite risky to insist a priori that dynamical chaos reproduces all properties of ran- 
dom fiows required for analytical approaches. I t is even less obv i ous that interstellar 
turbulence driven by supernova explosions (e.g. iKorpi et al.lll999l : ICressel et aLluOOSI ) 
provides a truly random velocity field. On the other hand, dynamo models based on 
steady fiows with stochastic flow lines such as ABC fiows excite magnetic fields which 
look rather different from the field discussed for turbulent dynamos. Zeldovich (see 
IZeldovich et al.ll 19831 ) suggested that dynamo action in nonstationary fiows (say, when 
parameters A, B and C in the ABC fiow fiuctuate in time) becomes much more similar 
to the dynamo action i n random fiows than dynamos in stationary fiows. Recent work 
of lKleeorin et al.l ( 20091 ) supports this idea and demonstrates that in a fiuctuating ABC 



flow, a large-scale magnetic field can indeed grow in a way similar to what is supposed 
to grow in a random mirror-asymmetric (helical) turbulent fiow. 



4 Conclusions 

In this review we have discussed our current understanding of both large-scale and 
small-scale dynamos that are relevant in astrophysics. In particular, we have illus- 
trated differences and similarities between them and have compared them in terms of 
the energy spectra with the corresponding cases at low magnetic Prandtl numbers. 
We have also briefly highlighted the resistively slow saturation phenomenon as well 
as catastrophic quenching of helicity-driven large-scale dynamos. Finally, we have dis- 
cussed new issues in connection wit h small-scale dynamos and their intermittency. 

Relative to earlier reviews (e.g. iBrandenburg and Subramanianll2005al ) there have 
been some unexpected advances regarding the nature of magnetic helicity fluxes and 
direct observational evidence for magnetic helicity in the solar wind. Another com- 
pletely unexpected development concerns the numerical detection of the negative ef- 
fective magnetic pressure instability in simulations of strongly stratified turbulence. 
It is expected that these developments will contribute to an improved understand- 
ing of the magnetic field generation in astrophysical bodies. On the technical side, 
IBrandenburg and Subramanian (2005a) discussed just the basics of the basics of the 
test-field method, but now this technique has provided significant insights into issues 
such as non-locality in space and time, as well as the nonlinear quenching of dynamo 
coefficients. 

In this review we have barely touched upon applications to actual astrophysical 
bodies. In fact, a lot of progress has been made by trying to model the Sun. Direct nu- 
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meric al simulations of convection in spherical shells has shown sisns of cy c lic larRe-scale 
fields JBrown et al.ll20f0l . [20nl : iKapvla et~all l20f ol : Ichizaru et al.ll2O10l : iRacine et all 



2011 ), but only for systems that are rotating at least 3 times faster than the Sun; see 



Brown et alj 1 201 j ) for simulations with otherwise realistic solar parameters. Similar 



results have also been obtained for local simulations of the galactic dyna mo, which only 



appe ars to be excited when the rotation speed is artificially enhanced ( Gressel et alJ 
[200«). This might well indicate that one is on the right track, but that the turbulence 
present in the system is exerting too much effective diffusion owing to it being domi- 
nated by rather large eddies, ft is conceivable that in reality the turbulent eddies would 
be smaller, lowering thereby the effective turbulent diffusion, which can at the moment 
(with the larger eddies) only be emulated by adopting faster rotation. Similar results 
have recently also been seen in models of the negative effective magnetic pressure in- 
stability, where direct numerical simulations showed the development of the instability 
only when there were enough turbulent eddies in t he domain and thereby th e turbu- 
lent diffusivity sufficie ntly small on that s c ale; s ee iBrandenburg et al.l (120111 ) and, in 
particular. Figure 17 of lBrandenburg et al.l (2013). Thus, much has been learned about 



turbulent dynamos and their relevance for astrophysical systems but as usual, much 
remains to be elucidated. 
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